NPS55-86-026 

NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




TECHNICAL 



A PARTIAL EVALUATION OF THE INTEGRATED 
TACTICAL DECISION AID (ITDA) SYSTEM 



JAMES No EAGLE 



DECEMBER 1986 



Approved for public release; distribution unlimited. 
Prepared for: 

Naval Air Development Center 
Warminster, PA 18974-5000 

FedDocs 
D 208.14/2 
NPS-55-86-026 




V ) 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



i ; " rj_ a i 



Rear Admiral R. C. Austin D. A. Schrady 

Superintendent Provost 



The work reported herein was supported in part by the Naval Air 
Development Center. 

Reproduction of all or part of this report is authorized. 

This report was prepared by: 



Unclassified 




DUDLEY KNOX LIBRARY 

NAVAL PQS TGRAnjM R J_, . _ 


SECURITY C-ASS t cat , 




REPORT 


DOCUMENTATION PAGT NrEREV CA 93943-5lo7 


la REPORT SECURITY ClASSIF.CA 

UNCLASSIFIED 


ib RESTR C vE MARK NOS 


2a SECURITY C.ASS F Ca 7 ON A^> 


S* 3 .B' 'ON A. A A3. Of RE»OR' 

Approved for public release; distribution 


2b DECLASSIFICATION DOWNGRADING SCHEDULE 


uni i mi ted. 




4 PERFORMING ORGANIZATION REPORT NUMBER(S) 

NPS55-86-026 



5 MONITORING ORGANIZATION RE POR r NUMBER'S, 



6a NAME OF PERFORMING ORGANIZATION 

Naval Postgraduate School 



6o OFFiCE SYMBOL 
(If applicable) 

Code 55 



7 a NAME Of MON TOPING ORGANIZATION 

Naval Air Development Center 



6c. ADDRESS (City, State, and ZIP Code) 



7b ADDRESS (City. State, and ZIP Code) 



Monterey, CA 93943-5000 



Warminster, PA 18974-5000 



8a. NAME OF FUNDING SPONSORING 
ORGANIZATION 


8b OFFiCE SYMBOL 
(If applicable) 


9 PROCUREMENT 


.NSTRUMENT IDENT F CAT ON NUMBER 


8c. ADDRESS (City, State, and ZIP Code) 




10 SOURCE OF FUNDING NUMBERS 






PROGRAM 
ELEMENT NO 


PROJECT 

NO 


TASK 

NO 


WORK oNlT 
ACCESS ON NO 



11. TITLE (Include Security Classification) 



A PARTIAL EVALUATION OF THE INTEGRATED TACTICAL DECISION AID ( ITDA) SYSTEM 



12. PERSONAL AUTHOR(S) 

Eagle, James N. 



13a. TYPE OF REPORT 


13b time COVERED 


14 DATE OF REPORT {Year, Month. Day) 


15 PAGE COu\* 


Technical 


FROM TO 


December 1986 


20 



16. SUPPLEMENTARY NOTATION 



17 


COSATI CODES 


18 SUBJECT TERMS ( Continue on reverse if necessary and identify by block number ) 


FIELD 


GROUP 


SUB-GROUP 


ASW, model, ITDA, search 



19 ABSTRACT (Continue on reverse if necessary and identify by block number) 



The mathematical accuracy and modelling reasonableness of the barrier and area 
search modules in the Integrated Tactical Decision Aid (ITDA) system were examined. 
Several problems were identified and solutions suggested. The use of ITDA was 
monitored during a 16-day sea period on the USS Carl Vinson (CVN70). 



20 DISTRIBUTION / A VAIL A 81 LI TY OF ABSTRACT 
Q UNCLASSIFIED/UNLIMTED □ SAME AS RPT 


□ DTiC USERS 


21 ABS t RACT SECURITY CLASSiFICA' 

Unclassified 


'ON 


22a NAME OF RESPONSIBLE INDIVIDUAL 

James N„ Eagle 


22b TELEPHONE (Include Area Code) 

(408)646-2654 


22c OFF i_E , v V30 

Code 55Er 



DD FORM 1473, 84 mar 



83 APR ed t.on may be used until exhausted 
ah otner editions are obsolete 



SECURITY CLASSiF’CA^ ON QP 



1 A ' - 

986 606 243 



O U.S Government P' 



fl Partial Li'aluation of the Integrated lactical Decision Hid liiUfli System 

by 



James N. Eagle 

Department of Operations Research 
Naval Postgraduate School 
Monterey. Cfl 93943 

Introduction 

The purpose of this technical report is tuuof old: 

f to critically examine the ITDfl FiSlU barrier and area search 
models for mathematical accuracy and modelling reasonableness, and 
2. to report on the use of ITDfl programs on 
USS CARL VINSON (CVN 701 during a id day period of high intensity, 
exercise operations (RIMPflC 86). 

Available for examination mere ITDfl Ver. 1 (including the BASIC 
language source code), the more recent ITDfl Ver. 2 (without source 
code), and the JOTS II Technical Reference Manual (Ref. I). Reference is 
made in this report to Version 1 source code to explain the functioning of 
the Version 2 program. This mas required because the Version 2 source 
code mas unavailable. Fortunately, the barrier and area search modules 
appear functionally very similar, if not identical, in both ITDfl versions. 

So the lack of Version 2 source code is not judged to be especially 
critical to this examination. 
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fiSLU Barrier Search Model 



Brief Description of Model . The ITDR RSUJ barrier search model is 
accurately described in Fief. 1. Its description here is primarily to 
establish notation. 

The barrier search model estimates the probability of detection of a 
target submarine penetrating a back-and-forth RSUU barrier (Figure 1.). 
The barrier penetration is considered in "target-stationary relative 
space". That is, the target is assumed to be stationary, and all relative 
speed for the encounter is provided by the searcher. The searcher's 
relative speed component across the barrier front is 

V s costa), 

where V s is searcher speed, and a is the barrier advancing angle (when 
,.v>0] or retiring angle (when a<0). The component of relative speed 
perpendicular to the barrier is 

V T + V§ sin(cd, 

where V T is target speed. 



Target window, UJ 



i 




Figure 1. Barrier Geometry 



For barrier penetration points uniformly distributed across a target 
window (UJ}.. the probability of the target coming within range R of the 
searcher (P D ) is the ratio of the area "covered" by the searcher in 
relative space to the total area that could be occupied by the target. In 
Figure 2. the covered area is shaded. In the ITDfl barrier search model, 
this area is obtained by numerical integration. Allowance is made for 
sensors which sweep less than 369° by reducing the detection range R. 
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Figure 2. Calculation of P 0 



Discussion . Although the modelling concepts used are good., there are 
some minor problems in the implementation. 

I. The calculated F' D is too large for target luindotu sizes (UJ) less than 
2(B/2 - RK). This problem results because the limits of the numerical- 
integration are, in some instances, improperly set. Specifically, in line 
51970 of the ITDA Ver. 1 BASIC code, the upper limit of integration is 
B/2 - RK. It should be min (B/2 - RK, IJJ/2). With the original limits, the 
integration can be performed outside the target window, W. This 
increases the calculated F' D above the proper value. To illustrate with 
an example, let 
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R = detection r ange = 10 nrn 
V s = searcher speed = 19 knots 
V' T = target speed = 19 knots 
B = barrier midth = 89 rim 
IJJ = target window = 19 nrn 
oc = advancing/retiring angle = 9 

K = kinematic enhancement factor = (1 + I. V 3/ V j j ^ 1 1 ~ = 1.414 
The ITDfi model gives a probability of detection of I. However, for any 
UJdB - 2RK), the correct answer is 2RK/B - .35. 

2. P D is improperly calculated when either V s cos(ct) or V T - V s sin(a) 
is negative. The cause of this problem is that the integrand in the 
numerical integration in line 51169 of the BASIC code is allowed to 
become negative whenever 

Ksqr = V s costa) / (V T + V s sinta)) < 9. 

Adding negative terms in the numerical integration results in the 
calculated P D being too small. To illustrate, continue the previous 
example with IU = 199 nrn and a = -19°. Then P D is correctly calculated as 
.31. Setting a = - 179° should give the same result, but yields .29. 

The failure becomes particularly spectacular when the denominator 
of Ksqr is exactly 9. For example, setting 
R = 19 nrn 
V' s = 29 knots 
V t = 19 knots 
B = 89 nrn 
ill = 199 nrn 
a = -39° 

results in a P D of 1.42 x I9 b 
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R simple fix for this problem is only allow <>. such that 
~90° ■■ ct 90° arid Vj + vg sin f oc } s 0 . 

The other cases are either redundant or result in the searcher moving 
"south" faster than the target. 

3. Sweeping sensors are modelled in an unusual manner The user is 
given the option of calculating P 0 under the assumption that the 
detecting sensor sweeps ±8° either side of a center bearing. The 
effective detection range then becomes Rsin(minl6,90°)). The problem 
with this approach is that it assumes that the center bearing is along the 
relative track of the searcher. Rn more reasonable assumption is that 
the center bearing be along the searcher's absolute (i.e... geographic) 
track. See Figure 3. 




Figure 3. Orientation of Swept Hrea 



These two assumptions lead to significantly different calculated 
probabilities of detection. For example, in the limit as 8 and ct decrease 
to 0 , the ITDR model gives a probability of detection of 0 . Whereas 
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assuming that the sensor is looking directly ahead along a fixed line of 
bearing gives a probability of detection of mini I, Fl/UJ). Note that in this 
case, there is no increase in Pq as the searcher increases speed li e , no 
kinematic enhancement). The probability of detection is simply that 
fraction of the possible barrier penetration points that are covered by 
the finite detection "ray" of the sensor. For nonzero values of 8 and ot, 
similar closed form expressions can probably be derived, although the 
geometry becomes more complex and performing the numerical 
integration may be preferred. 

4. This is a comment, not a problem. For c*=0° and a 369° sensor, it is 
possible to approximate the ratio of the shaded area to the total area in 
Figure 2. without numerical integration. True, the HP9929 is not being 
computationally stressed by the integration, but there is something to 
be said for the simplicity of closed form solutions. Using the notation 
introduced, the probability of detection given a barrier penetration is 
approximately 



min (2RK/B, 1J for LU $ B — 2RK 



Po = 



1 - (1/UJB)(max (,5(B + UJJ - RK, 0jj 2 for B-2RK < UJ < B + 2R 

((B + 2RJ/UJ) (1 - (max (B - RlK-1), 0)J 2 /(B 2 +2RBJ) for IJJ > B + 2R 



where K = (I + (V 3 /V T ) 2 ) ,/2 as before. To obtain this expression for F' 0 , 
the circular boundary of the shaded area in Figure 2. is replaced with a 
straight line The resulting region has an area which is greater than or 
equal to the area of the original shaded region So this calculation gives 
an upper bound for P 0 . Figure 4. is a comparison plot of P D calculated by 
the ITDfi model and the approximation above The ITDfi model values 



greater than .35 are in error and result from improper limits of 
integration used when calculating P D for small IU. 




Figure 4. Comparison of ITDR Model Results and Approximation 



5. This is a second comment. One advantage of using numerical 
integration for evaluating P D is that it can accommodate a nonuniform 
distribution of target penetration points. This advantage is not 
exploited in the current barrier search program. For example, if the 
target knows that the barrier is in place, then an attempt might be made 
to penetrate near the one of the ends of the barrier to minimize F' 0 . If 
numerical integration is used to evaluate P D , then it is a simple matter to 
weight the possible barrier penetration points according to any specified 
probability distribution. Specifically, if the barrier penetration point x 
has a probability density function g(x) and P D (xl is the probability of 
detection given a barrier penetration at x, then 



U 



p d = J‘ p d (x) gtxi dx 

HI I ;< 

s.t q-'x»0 

6 R third comment. R possible problem ujith this advancing/retiring 
barrier methodology is that the detection probabilities are conditional on 
the target's penetration of the barrier Changing the angle of searcher's 
advance or retreat ijuill affect both the probability of detection given 
barrier penetration and the probability of barrier penetration. The 
second effect is not accounted for in this model. For example, a 
retreating barrier can dramatically increase the probability of detection 
given a barrier penetration, but it also reduces the probability of such a 
penetration during any specified time interval. In the limit as the 
searcher speed perpendicular to the barrier increases to the target 
speed, the predicted P D increases to 1. However, the target can not now 
overtake the retreating barrier, so the probability of a barrier 
penetration occurring at all becomes zero. This line of argument 
suggests that the calculated probabilities of detection may be too large 
for retiring barriers and too small for advancing barriers. 

flSUJ Area Search Model 



Brief Description of Model , fls described in Ref. 1, probability of 
detection for an area search is given by 



iij 



P D = I - expl-ZTi V R t/A), 

inhere 

R = detection range 
fl = search area size 
t = time alloiued for search 
V's = searcher speed 
V T = target speed, and 
V R = mean relative speed 

2p 

= 11/2 rr I j (V.-/ + V T 2 + 2V s V t cos 4>) 1/2 d<t> 
u 



(2 



Discussion 



1. The mean relative speed calculation can be made more efficient. 
In the original BASIC code and in Ref. 1, the calculation of V R is as shown 
in Equation (2). A preferred expression for numerical evaluation is 

Tl 



(1/ JT ) | 



- 2V s V t cos 4') 1/2 d$ . 



(3) 



Equation (3) makes use of the symmetry of the cosine function to reduce 
the extent of the numerical integration, and has the correct sign for the 
2 V s V t cos'4 term in the Lam of Cosines. These problems are to some 
extent cosmetic, as cosine symmetry causes both (2) and (31 to give the 
same numerical results. 

It is noted that V R can be approximated mith the following closed 
form expression: 

V R /V S ~ .36 max(V T /V s , 1) + .64 (1 + (V T /'V S ) 2 ) 1/2 (4) 

Figure 5 contains plots of the estimated V R /V S versus V T /V S and a 
multiplicative correction factor to convert the estimated values to true 
values. The estimate is always within 1.1% of the true value. 
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SPEED RATIO (TGT. SPD./SEARCHER SPD.) 




Figure 5 Estimated Relative Speed Ratio vs Speed Ratio 



2. R Monte Carlo area search simulation mas conducted to test the 
accuracy of Equation 11). The search geometry is shown in Figure 6 The 
searcher starts in the lower Left corner of a IRQ rirn x IQ 9 nrn search area 
arid conducts a systematic search at speed V s with a cookie-cutter 
sensor having a range of IQ nrn. The target's starting position is 
uniformly distributed over the search area. The target conducts a 
random tour at speed V T with a mean time between course changes of 1 
hour. That is, the target track consists of connected Line segments, 
each of which has a course selected from an independent, uniform 
distribution between 0° and 360°. find the time on each leg is selected 
from an independent, exponential distribution with mean 1 hour. When 
the target encounters a search area boundary a perfect reflection is 
performed. 



1 00 nm 




Figure 6. Simulation Geometry 



Seven simulations mere conducted, each mith V s and V T selected so that 
the mean relative speed, as calcuated from Equation 13), mas 15 knots. 
The results are plotted in Figure 7, together mith P 0 calculated from 
Equation I'D. The fit betmeen (1) and the simulation mas reasonable when 
Vc and V T mere nearly equal. Otherwise the fit was not especially 
impressive. 

The simulation results were sensitive to the mean time between 
target course changes and the searcher track used (neither of which are 
represented in Equation (!)]. In particular, the simulation showed that a 
search track starting near the center of the search area gives better 
search performance, especially for small values of V s . 

3. The need for an area search simulation. For this area search 
example. Equation (D does not model a systematic search for a randomly 
moving target very accurately, find unfortunately the search literature 
has little else to offer in the way of closed form expressions for P D 
when the target track is random and the searcher conducts a systematic 
search (e.g., parallel sweeps, spiral in, or spiral out). This suggests that 
some sort of numerical solution, probably a Monte Carlo simulation, may 
be required to assess the effectiveness of such an area search. This 
simulation could be very general. It could allow, for example: 

1. different areas for the target and searcher, 

2. a customized target motion model (e.g., patrol or transit) 

3. any desired search track, 

4. multiple searchers and targets, 

5. counterdetection by the target. 



PROBABILITY OF DETECTON 



Mean Relative Speed is 1 5kts 
Mean Rate of Course Changes is 1/hr 



o 




Figure 7. Probability of Detection by Simulation 
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Output could include the searcher trad Is) and graph of probability of 
detection by time t. With a cursor or light pen input of area boundaries 
and track points, problem setup time cuould be reduced This simulation 
might be of particular use for battle group flSW planning, since 
frequently several FiSlU platforms search for the same target 
submarine 



fit Sea Observations 



Back g round During the period 2-18 June 1986, the author, as a 
commander in the Naval Reserve, mas attached to the staff of 
Commander Carrier Group Three, fit that time, the staff mas embarked 
on USS CARL VINSON (CVN 70) and operating in support of the exercise 
RIMPRC 86. 

RIMPfiC 86 mas a two-carrier battle group, multinational exercise 
involving approximately 60 ships and 30,000 at-sea personnel The 
purpose mas to demonstrate the fleet's ability to operate effectively 
against a realistic air, surface, and subsurface threat. During this 
exercise, use of the ITDH system by the carrier group staff and CARL 
VINSON ship's company mas observed 

Exercise Scenario. The general scenario mas that a coup had toppled 
the pro-West government of White Country (represented for purposes 
of the exercise by the Hawaiian Islands) The US naval response mas to 
send two carrier battle groups (Ranger and Carl Vinson) to White 
Country as a shorn of force. These battle groups mere opposed by 
some 60 ships simulating Soviet platforms. As the exercise progressed. 
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the battle groups approached iJJhite Country close enough to launch air 
strikes. Opposing surface action groups and submarines engaged as the 
battle groups came into range. 

JOTS and DOTS, The ITDR system on CRRL VINSON is called VINSON 
DOTS, which appears to be a one-of-a-kind installation, funded by 
ORRF'fi, through NOSC, San Diego. There was no formal documentation 
for this system aboard CRRL VINSON, but the overall program structure 
and functional modules seem very similar to the older JOTS. In fact, 
much of the BRSIC code has been carried over directly from JOTS (as 
indicated by documentation comments in the BRSIC code). The principle 
DOTS enhancements included the organization of the contact database 
and the implementation of automated input of LINK 11 data. 

The controlling terminal for VINSON DOTS is an Hewlett-Packard 9020 
located in the UJar Room of the Flag Spaces. Feeding this station is 
another HP 9020 in the Carrier Intelligence Center (CVIC). This unit 
processes and correlates intelligence information generated off the 
ship. The correlator-tracker used is POST (Prototype Ocean Surveillance 
Terminal). Also feeding into the Flag HP 9020 is real-time NTDS (Link 11) 
data. The Flag terminal displays and operates on available contact 
data and also sends this data to HP 9020 machines in the RSLU Module 
(RSLUMQO) and Combat Direction Center (CDC). Another HP 9020 is located 
in the Tactical Operations Plot (TOP) which runs JOTS III and is currently 
not in the VINSON DOTS network (and so does not have access to real- 
time NTDS information). 



The VINSON DOTS network accepts data only from the Flag terminal 
and these data currently come from three sources: II manual input, 2J 
NTOS input, and 3) POST input. 

The NTDS input is automated and requires no human filtering 
Automated input is also available from POST, but due to the 
classification of some POST data, this feature is not presently 
implemented. To remove the more sensitive data, POST contact 
information must be manually sanitized and then entered into the 
VINSON DOTS database from the Flag terminal. Automatic sanitization 
and entry of POST data is a planned enhancement. 

Even though contact data to be shared over the network must be 
entered from the Flag 9Q20, any other terminal may be used to enter or 
modify contact information to be used exclusively on that terminal. 

This allows users to customize their databases. 

Use of DOTS Aboard CVN 78. Of the many tactical decision aids that DOTS 
provides, the principle observed uses during RIMPAC 86 were. 

1. Maintaining an contact history database for the exercise, 
containing all surface, air, and subsurface contacts. 

For each detection recorded, some or all the following data were 
saved: time, location, target course and speed, detecting sensor and 
detection parameters (dimensions and orientation of elliptical spas or 
bearing spread and inclusive ranges for bearing boxes). In addition, 
■JOTS III (but not DOTS) provides a comment field which was used 
extensively by the TOP watchstanders (responsible for maintaining the 
surface plot). 

2. Plotting current or past contact positions. 



There were two principle ways this capability was used: 

a. Displaying the current position of some group of contacts le g.. 
Blue surface ships), or 

b Displaying the contact history (i.e., track) of a few contacts. 

Color and NTDS symbology were effectively used to identify the 
platform type and threat status (hostile, friendly, or unknown). 

The processing and display of Link 11 data transfer was fast enough 
to follow the air war. This was demonstrated during several air raids 
on the battle group where actual aircraft were used, (fl-3's simulated 
bombers with A-4's and fi-7's flying as air-launched cruise missiles.) 



- e.g., if BATES was underway from San Diego 9890 on 2 June, could 
she be in the Hawaiian Operation area by 9809 on 10 June with a 15 knot 
speed of advance? 

r with a 18 knot speed of advance going through specified points, 
where will the contact be in 24 hours? 

4. Plotting spas, patrol, areas, 4-lUhiskey grid (i.e., screen) 
assignments 

5. Preparing "executive summaries" of the contact situation. 



No direct use of the ARID, ASUI or ELU decision aids was observed, 
although several of the Carrier Group Three staff members were very 
knowledgeable in their use. 



0 b s e i ' v 9 d P ro b 1 e rn s . 

1. Resolution on the HP 9828 screen was insufficient to display many 
contacts. (The POST display, using an HP 98280 machine, was much 
better.) 

2. TOP watchstanders used JOTS III rather than DOTS because JOTS 
III is more complete. 

3. The HP Think Jet printer generally only provided marginally 
acceptable hard copy. The installed thermal printer on some of the 
HP 9828 machines provided a finer resolution output. 



References 



(1) JOTS II Technical Reference Manual . Commander Naval Air Force, U S 
Atlantic Fleet, May 1984, CONFIDENTIAL. 



19 



DISTRIBUTION LIST 



NO. OF COPIES 



Defense Technical Information Center 2 

Cameron Station 
Alexandria, VA 22314 

Li brary 
Code 0142 

Naval Postgraduate School 
Monterey, CA 93943 

Research Administration 
Code 012 

Naval Postgraduate School 
Monterey, CA 93943 

Li brary 
Code 55 

Naval Postgraduate School 
Monterey, CA 93943 

Center for Naval Analyses 
2000 Beauregard Street 
Alexandria, VA 22311 

Operations Research Center, Room E40-164 
Massachusetts Institute of Technology 
Attn: R. C. Larson and J. F. Shapiro 

Cambridge, MA 02139 

Koh Peng Hong 1 

Mi ni stry of Defence 
Blk A, Stockport Road 
Singapore 0511 

Naval Air Development Center 1 

Attn: Dr. M. Leonardo 

Warminster, PA 18974 

Naval Air Development Center 1 

Attn: LCDR R. Hudson 

Warminster, PA 18974 

Navy Tactical Support Activity 1 

Attn: Mr. A. Letow, Director 

P.0. Box 1042 

Silver Spring, MD 20910 



2 



1 



1 



1 



20 



mini FY KNOX LIBRARY 




3 2768 00302445 6 



